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The ghost propagator and the Coulomb potential are evaluated in Coulomb gauge on the lattice, 
using an improved gauge fixing scheme which includes the residual symmetry. This setting has 
been shown to be essential in order to explain the scaling violations in the instantaneous gluon 
propagator We find that both the ghost propagator and the Coulomb potential are insensitive 
to the Gribov problem or the details of the residual gauge fixing, even if the Coulomb potential 
is evaluated from the Ao-propagator instead of the Coulomb kernel. In particular, no signs of 
scaling violations could be found in either quantity, at least to well below the numerical accuracy 
where these violations were visible for the gluon propagator. The Coulomb potential from the 
Ao-propagator is shown to be in qualitative agreement with the (formally equivalent) expression 
evaluated from the Coulomb kernel. 



8th Conference Quark Confinement and the Hadron Spectrum 
September 1-6, 2008 
Mainz. Germany 



* Supported by DFG under contract Re-856/6-1,2 
^Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Coulomb gauge ghost propagator and the Coulomb potential 



Markus Quandt 



1. Introduction 



Yang-Mills theory in the Coulomb gauge has recently drawn a renewed attention, both in the 
continuum [[I], |, ||] and on the lattice [|, H, ^ H). This is mainly due to the fact that GauB' 
law can be resolved explicitly in this gauge, which allows for a neat Hamiltonian formulation with 
the transversal part of the remaining vector potential as the only physical degree of freedom. 
Much of the intuition and techniques from ordinary quantum mechanics can thus be carried over 
to the YM case. In particular, recent variational approaches in the Schrodinger picture, based 
on the notion of a weakly interacting constitutent gluon and the Gribov-Zwanziger confinement 
scenario [12], proved to be very successfuU [^; similar calculations are presently carried out in the 
renormalisation flow approach. 

All these continuum formulations, in one way or the other, give rise to relations between low- 
order Green functions of the constituent gluon A and the Faddeev-Popov ghosts. It is therefore 
important to obtain non-perturbative information on such correlators from the lattice. Careful stud- 
ies of the equal-times gluon propagator, for instance, reveal strong scaling violations and a UV 
behaviour at odds with simple dimensional arguments [§, 0, ^. These surprising results reflect the 
renormalisation problems for instantaneous correlators in the continuum. One possible explana- 
tion of the lattice findings is based on the idea that the residual gauge freedom left over by the 
Coulomb condition must be fixed in such a way that it resembles the Hamiltonian formulation as 
closely as possible. ' A careful study of the energy dependence of the gluon propagator then allows 
to manipulate the data such that perfect scaling is observed even on finite lattices. 

For the confinement scenario layed out by Gribov and Zwanziger[12], the more important cor- 
relators are, of course, the ghost propagator and, in particular, the Coulomb potential. Furthermore, 
the ghost form factor has been shown to represent the inverse of the colour dielectric function of the 
Yang-Mills vaccum [13|, and is therefore of direct physical relevance. Initial studies of the ghost 
and Coulomb propagator for the gauge group G = SU (2) with simple Coulomb and no residual 
gauge fixing ^ found no scaling violations at low momenta, but had inconclusive results about the 
Coulomb string tension in the deep infrared. Moreover, these results were partially at odds with 
more careful SU (3) studies using a residual gauge fixing different from ours [11], which featured 
a peculiar saddle-like behaviour in the Coulomb potential at low momenta. In the present talk, I 
will report about recent SU (2) calculations of ghost form factors and the Coulomb potential, using 
exactly the same gauge fixing techniques which proved essential for the resolution of the scaling 
violations in the gluon propagator. 



2. Gauge Fixing 

Our gauge fixing procedure employs both simulated annealing and the microcanonical flip 
procedure layed out in [^] as a preconditioning with subsequent (over)relaxation to complete the 
gauge fixing within machine precision. To reduce the Gribov noise and bring the lattice configs 
closer to the fundamental modular region, we perform up to 40 restarts with random gauge transfor- 
mations as starting points, and take the copy with the best minimum of the gauge fixing functional. 
While this procedure proved to be important for the correct extraction of the gluon propagator in 

'For the first-order formalism in tlie continuum, renormalisability has been proven algebraically IlKJ]. 
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Figure 1: Left panel: Energy dependence of the Ao-propagator Doo(p,/?o) after improved Coulomb and 
resiudal gauge fixing, for various spatial momenta |p|. Right panel: The equal-times Ao-propagator 
Doo{p,t = 0) as a function of the spatial momentum |p|. 



the deep infrared [0], the ghost correlators exhibit a much weeker dependence on the quality of 
gauge fixing. This can be clearly seen in the left panel of fig. § The value of the ghost propagator 
at the lowest diagonal momentum p = (1,1,1) is only very slightly suppressed as the number n of 
Gribov restarts is increased, and the optimum is already reached for n as low as n ^ 2.3. All this 
is in constrast to the corresponding findings for the gluon propagator, where a 20% effect was seen 
that required up to « = 40 for saturation. 

The second important ingredient is the residual gauge fixing. To make contact with the Hamil- 
tonian approach in Weyl gauge, we would like to put the spatial average u{t) = L^^^x^o(?,x) 
to unity. However, periodic boundary conditions only allow us to make u{t) time-independent, 
u{t) = Uq = const. In the infinite volume limit (and in praxis also for L > 32), Uq approaches 
unity. Although this only enforces ^o^^'o = on the spatial average, the Ao-propagator is, within 
statistical errors, independent of energy (see left panel of fig. |l]). In the right panel of fig. |l|, we thus 
plot only the instantaneous Ao-propagator which is strongly enhanced in the infrared. This result 
will be related to the Coulomb potential below. 



3. Results 

The right panel of figure ^ shows our results for the ghost propagator and its form factor, 

G{p) = (c-(-p)c(p) ) = £ e'P^ ( M(x,0)-' ) ^ ^ (3.1) 

X P 

where M = (— VD) is the Faddeev-Popov operator and the ghost form factor d{p) measures the 
deviation from the perturbative result. The form factor is infrared enhanced, which agrees with the 
horizon condition d^^(0) = necessary in the Zwanziger confinement criterion [12]. Our infrared 
exponent K ^ 0.22 for the divergence d{p) ^1/ {p^)^ is slightly smaller than the one obtained with 
naive gauge fixing [Q], but agrees well with recent improved studies in SU (3) [jlT]]. 

Even more directly related to the confinement problem is the so-called Coulomb potential Vc, 
i.e. the response of the gluon vacuum to static colour charges. Since the constituent gluon A and 
its wave functional are gauge-dependent, Vc is not directly the physical potential between static 
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Figure 2: Left panel: The ghost propagator at the lowest diagonal momentum p = (1, 1, 1) as a function of 
the number of Gribov copies considered in the Coulomb gauge fixing. (Note the scale on the y-axis.) Right 
panel: The ghost form factor d{p) as a function of the spatial momentum |p|. 



quarks (as extracted from Wilson loops or Polyakov lines), but an upper bound, Vc{r) > jV{r) 



This implies that there is no confinement without Coulomb confinement []12|], but a linear Coulomb 
potential may persist even in the deconfined phase. 

Formally, Vc (r) can be computed in one of two equivalent ways. 



V,(|x-y|) = (Ao(r,x)Ao(r,y))=g2((M-i.A.M-i) ). 



(3.2) 



The formal equivalence of these two expressions can be shown in the first order formalism upon 



explicitly resolving GauB' law |1C, 14]. This leaves possible renormalisation issues aside and 
the lessons learned from the scaling violations in the gluon propagator indicate that some caution 
is required when connecting bare instantaneous correlators. Of course, the Ap-propagator is nu- 
merically much simpler than the complicated Coulomb kernel involving two inversions of the FP 
operator. 

The strong Ward identities in Coulomb gauge [ p^ imply that the special combination in mo- 
mentum space 

?Hp) (3-3) 



P'K:(P) 



is a renormalisation group invariant which can be taken as a definition of the running coupling con- 
stant. Simulations with different j3 should thus fall on top of each other without further multiplica- 
tive renormalisation. We have tested this conjecture for numerous values of j3 on relatively small 
16^ lattices. (On 32'* lattices, we have only been able to complete the analysis of the complicated 
Coulomb kernel for a single value of j3). The jS-invariance was much better for the Ao-correlator, 
while Vc constructed from the Coulomb kernel still showed noticable scaling violations. At present, 
it is not known whether these deviations are pure numerical or finite volume effects, or if they have 



any more significant meaning. (Similar observations were made in ref. flllQ). Simulations with 
improved statistics on larger lattices have to be conducted to resolve this issue. 

Finally, the most direct approach to the confinement issue is given by the expression 



p'V,(p). 



(3.4) 
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Figure 3: Left panel: The combination Vc (|p|) with the Coulomb potential Vc extracted from the Aq- 
propagator Doo(p,f = 0). Right panel: The same quantity, with Vc extracted from the Coulomb kernel. 



From the Fourier transformation of a linear potential, Vc{r) = Oc r, it is readily seen that 

pV,(p)^87rcT,, IpHO. 

The Coulomb string tension Oc is an upper bound for the real string tension a extracted from 
Wilson loops. Previous and current lattice studies are inconclusive as to whether = a, since the 
approach to |p| — > is not as uniform as expected: Early simulations without improved/residual 



gauge fixing saw a slight but noticeable rise in the quantity (3.4) below |p| ~ 1 GeV, which seemed 
compatible with Ocjo anywhere in the range 1 . . .3. More recent computation for the gauge group 
G = SU (3) prefer a value ac/cJ « 1.6, but the extrapolation to zero momentum is again uncertain 
due to a peculiar "bump" in the quantity ( |3.4[ ) at momenta between 0.1 ... 1 GeV. 

Our result on aV = 32^ lattice in figure ^ using all improved gauge fixing techniques give 
reliable results (for cylinder cut momenta) only down to |p| ~ 0.5 GeV. In this range, the results for 



(3.4) are compatible with Vc computed either from the Ao-propagator or from the Coulomb kernel. 
The latter result show a more pronounced plateau at the smallest momenta, which is reminiscent 
of the slight rise observed in However, the numerical data can equally well be fitted with 
a constant. (Vc from the Ao-propagator is compatible with the Coulomb kernel results within 
statistical errors). For both definitions of Vc, we do not see the "bump" reported for SU (3) in 
ref. [pl|]. While the approach to a constant seems promising, better statistics and larger lattices are 
required for a reliable extrapolation of Oc/ CJ. 



4. Conclusion 

The computation of ghost correlators and the Coulomb potential in G = SU (2) show qualita- 
tive agreement with continuum calculations in the variational approach |@, |3|]. The scaling viola- 
tions observed previously for the equal-times gluon propagator D{p) have no counter part in the 
ghost correlators studied here. In particular, the dependence on the Gribov noise and the details of 
the improved gauge fixing are negligable. Likewise, the residual gauge fixing, which is essential 
for the resolution of the scaling violations in D{p), seems to have little or no influence on the ghost 
propagator or Coulomb potential, even when the latter is extracted from the the Ao-propagator. 
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Our residual gauge fixing removes the energy dependence on Aq not only in the spatial aver- 
age, but effectively for arbitrary Ao-correlators. There is thus no issue with renormalisation and 
the results for the instantaneous Ao-correlators resemble the ones with unfixed residual symmetry. 
(Similar observations are made for the ghost propagator and the Coulomb potential as extracted 
from the Coulomb kernel.) It is therefore not surprising that our findings agree with other calcula- 
tions, even if these fixed the Coulomb gauge naively, or left the residual symmetry unfixed. 

The statistics in the deep infrared are not sufficient to make reUable quantitative extrapolations 
for the Coulomb string tension Oc, or the Coulomb form factor f{p) whose infrared behaviour is 
an important ingredient in the variational approaches [Q, ^. We intend to improve on this and 
accumulate data for 32^ lattices with various j8, and Ao-correlators on even larger lattices. These 
results will be published in a forthcoming paper. 
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